Universal contact for a Tonks-Girardeau gas at finite temperature 
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We determine the finite-temperature momentum distribution of a strongly interacting ID Bose 
gas in the Tonks-Girardeau (impenetrable-boson) limit under harmonic confinement, and explore 
its universal properties associated to the scale invariance of the model. We show that, at difference 
from the unitary Fermi gas in three dimensions, the weight of its large-momentum tails - given by 
the Tan's contact - increase with temperature, and calculate the high-temperature universal second 
contact coefficient using a virial expansion. 
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Introduction. Ultracold atomic gases are a unique sys- 
tem to probe universality aspects in physics. Due to their 
diluteness condition, the interactions are in most of the 
cases modeled by a single parameter, the s-wave scat- 
tering length, which describes the only type of collisions 
allowed at the ultralow temperatures typical of the ex- 
perimental conditions. Another remarkable aspect is that 
the interactions are tunable (eg by the mechanism of Fes- 
chbach resonances). When the interactions are brought 
to extremely large - positive or negative - values, pro- 
vided that the system remains stable against collapse or 
three-body recombination, there is no energy or length 
scale associated to the interactions. In this regime, char- 
acterized by spatial scale invariance, all the properties 
of the system are governed only by external parameters, 
such as the temperature, the density or the external con- 
finement, and the system thus displays universal features. 
For the case of a three-dimensional two-component Fermi 
gas in the so-called unitary regime (corresponding to in- 
finitely large scattering length) several universal aspects 
have been explored both theoretically [H-Q and experi- 
mentally [10l4l5| . 

We consider here a one- dimensional (ID) Bose gas, 
tuned to the strongly repulsive regime, known as 
impenetrable-bosons or Tonks-Girardeau (TG) limit. 
This regime has been experimentally achieved [Mil. 
From the theoretical point of view, a mapping onto a 
Fermi gas [22] yields the exact many-body wavefunction 
of the system, allowing to obtain with extreme accuracy 
information on the state of the system and on its quan- 
tum dynamics. In fact, it is its fermionic character which 
ensures the stability of the system in the experiment [2^] . 
In this work we explore some of the universal aspects of 
the ID Bose gas in the TG regime. In particular, we 
focus on the Tan's contact coefficient @, which corre- 
sponds to a two-body correlation function. Being associ- 
ated to the average interaction energy in the system, it 
directly reflects the universal aspects of the system once 
the interactions are tuned to very large values. The Tan's 



contact can be extracted from the wings of the momen- 
tum distribution [a, Ylm . Momentum distributions have 



been accurately measured in the experiment [25j. 

A very important experimental issue is the effect of the 
temperature. The various finite-temperature re gim es of 
the ID Bose gas have been previously identified [26l428| . 
However, the momentum distribution of a ID Bose gas at 
finite temperature has received relatively little attention 
till now. For a homogeneous Bose gas it is in principle 
accessible by thermal Bethe Ansatz calculations by ex- 
tension of Ref. [2^] . From the large-distance properties of 
the one-body density matrix, the Luttinger-liquid model 
predicts a Lorentzian shape at small momenta [30j. The 
momentum distribution for a TG gas was studied in a 
seminal work by Lenard [31| , who was mainly interested 
in the thermodynamic limit. The momentum distribu- 
tion of a TG gas under both box trap and harmonic con- 
finement has been numerically evaluated through lattice 
simulations for hard-core bosons at low filling [32]. We 
obtain here for the first time an expression for the tails 
of the momentum distribution at finite temperature and 
extract the universal contact coefficient. Our approach is 
valid for any external confinement, and is mostly analyt- 
ical. As a main result, we obtain that the weight of the 
universal tails of the momentum distribution increases 
with temperature. 

One-body density matrix from the thermal Bose-Fermi 
mapping. We consider Nb bosons of mass m confined by 
the harmonic potential V(x) = muj 2 x 2 /2. The particles 
interact via the contact potential v(x) — gS(x), and the 
Hamiltonian is given by 
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j<£ 



XI) 



(1) 



In this work we focus on the impenetrable limit g — > 00, 
where the effect of contact interactions can be replaced 
by a cusp condition on the many-body wavefunction. 
The repulsions are so strong that the many-body wave- 
function vanishes when two particles meet, ie ^(...Xj = 



2 



xe...) — for each pair {j,i}- The exact solution of the 
many-body Schroedinger equation H^N.a = EN,a^N,a 
for any TV-particle state individuated by the quantum 
numbers a — ...fjv} satisfying the above cusp condi- 
tion is readily obtained by the Bose-Fermi mapping [33 1 . 
This allows to obtain the bosonic wavefunction in terms 
of the one of an ideal Fermi gas in the same external 
potential and with the same quantum numbers, times a 
mapping function which ensures bosonic symmetry under 
particle exchange, 



^N, a (xi...x N ) = U 1 <j < t< N eigp.(xj-Xi)^ a (xi,x 2 ..,x N ). 

(2) 

Here, ^^ a (xi,x 2 ..,x N ) = det^ (a;*)] is the 

fermionic wavefunction constructed with the single parti- 
cle orbitals u v .(x). For harmonic oscillator confinement, 
rescaling the spatial coordinate in units of the harmonic 



oscillator length a/ lo = ^/h/moj, the single particle or- 
bitals read u v {x) — H v {x)e~ x ' 'I 2 /\/2 v v\/tt 1 / a with corre- 
sponding single-particle energies e Vj = Hu}(vj + 1/2). 

The above Bose-Fermi mapping allows to construct the 
thermal average of any observable. Observables that do 
not depend on the sign of the many-body wavefunctions 
are readily given by their fermionic counterparts 33], 
while substantial differences arc expected for those ob- 
servables which do depend on the sign of \1/ as it is the 
case of the one-body density matrix. Its expression at 
temperature T in the grand-canonical ensemble reads 

PiB{x,y) = ^ Pn,<xN / dx 2 , ...dx N 

N,a ^ 

x ^N,a(x, x 2 ..,x N )^* Na (y,x 2 , -,x N ). (3) 

Here I — (— oo,oo) is the spatial integration domain for 

jZ is the 

-/3(E NiC ,-ij.N) 



harmonic confinement, Pjv a 



-/3(E N , a -fiN) 



N.a ' 



thermal distribution function, Z = 
the partition function for the Tonks- Girardeau gas (or the 
mapped Fermi gas) with E N ^ a = J2j=i e v v ft = 1/ksT, 
and [i the chemical potential. 

The evaluation of Eq.([3]) as it stands appears as a 
formidable task, but we simplify it as in an early work by 
Lenard [31( (see the supplementary material for details). 
The final result reads 

PiB(x,y) = ^2 — ~ (sign^ _ v)) j / dx 2 ...d£Cj+i 



x det[p 1F (xi,xe)] iii=1 j + i, 



(4) 



where pif(x, y) = Y^j=i UjU Vj {x)ul. (y) is the fermionic 
one-body density matrix, f v = \/\e^ ev ~ IJ '' + 1] is the 
Fermi occupation factor of a single-particle energy level, 
and in the above determinant one has to take Xi — x for 
i = 1 and xi = y for t = 1. This illustrates how also for 
off-diagonal coherences of the Tonks-Girardeau gas it is 
possible to resort to the solution of a fermionic problem, 



thus more involved than the one needed for the density 
profile. 

Momentum distribution of a Tonks-Girardeau gas un- 
der harmonic confinement. Equation (j4|) is our start- 
ing point for the evaluation of the bosonic one-body 
density matrix, which, by Fourier transform, yields the 
momentum distribution of the gas according to n{k) = 
(l/2ir)fdxdye* k ^p 1B (x,y). 

First of all, we proceed to further simplifying Eq. (j4}. 
Expanding the determinant in Q and using the defini- 
tion of the fermionic one-body density matrix, the j-th 
term of the one-body density matrix according to the ex- 
pansion p 1B (x,y) = X\[-2sign(a; - y)] 3 Px B (x,y)/j\ is 
given by 

Pib( x 'V) = f^-fvj+i {- l ) Pu u 1 {x)u l/p{1) {y) 

vi..vj+i PeSj + i 



Y[ dx e u Ue ( 



(5) 



1=2 J y 

This can be finally casted onto the compact form 

i+i 



where 



-4,, 



fc=i 



Xx,y) = (-l) 1+k detB 



(6) 



(7) 



and detH UlVk (x,y) is the minor determinant of the ma- 
trix function B(x, y) taking out the row 1 and column 
k, whose full expression is given in the supplementary 
material. We have thus considerably simplified the com- 
plexity of the problem, by reducing a j-variable multi- 
dimensional integral to the combination of single- variable 
integrals, which moreover can be performed analytically 
in terms of special functions (see the supplementary ma- 
terial for details) . We are finally left to numerically evalu- 
ate only the sums over the single-particle levels in Eq. © . 

The momentum distribution of the TG gas at finite 
temperature is illustrated in FigfT] [42]. At increasing 
temperature we notice a decrease of the peak at small 
momenta^ (this was also observed in the numerical simu- 
lations [32]), while at small temperatures the wings of the 
distributions appear to be less affected by the thermal ex- 
citations. The high-momentum tails of the distribution 
are highlighted in the inset, where the momentum distri- 
bution is shown in double-logarithmic scale. As a main 
result, we find that the weight of the tails increases at in- 
creasing temperature. This is remarkably different from 
the case of a unitary three-dimensional Fermi gas. As we 
shall see below, this difference stems from the effect of 
reduced dimensionality. 

High-momentum tails at finite temperature The mo- 
mentum distribution of a ID Bose with contact inter- 
actions displays at high momenta a universal power-law 
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FIG. 1: (Color online) Momentum distribution of a Tonks- 
Girardeau gas (in units aho) as a function of wavevector (in 
units of a~^l) with N = 5 particles under harmonic confine- 
ment at increasing temperature, from top to bottom in the 
main peak k B T/fiw = 0.1,0.5,0.7,5,10,15. The last three 
curves at high temperature are obtained by the use of the 
classical fermionic one-body density matrix (|12|) into the se- 
ries expansion ((4]) up to j — 1. The inset shows the same 
curves in double- logarithmic scale. 



t£*F(z ) cos(7ra/2)r(a) for k 
the expression for the contact 

2 



oo we finally obtain 



C 



dRF(R). 



(11) 



At high temperatures, to estimate the high-momentum 
tails we use the classical limit for pip(x, y) 



Pif{x,v) - Aexp 



( * + y) -tanh-^ 



4 

(x-y) 2 



coth 



2 

jihw 



(12) 



with A = e^e"^"/ 2 /[^(l-e' 2 ^)] 1 ' 2 , to construct the 
j = 1 term in the sum (U) for piB(x,y). From its short- 
distance expansion, the classical limit of the two-body 
function reads 



F 



class 



A 2 Bhu 
( R ) - — coth—- exp 



-2i? 2 tanh 



(13) 



The large-temperature behaviour of the contact is then 
obtained from Eq. fTTj) and (fTBl as 



decay 35j, |36| n{k) — > C/k 4 for k — > oo, independent on 
interaction strength. The weight C is now known as the 
Tan's contact p], and is related to quantum average of 
the interaction energy Hint = 9 J2j<e ^( x j ~ x l)i 



C 



h 4 



(Hint) 



(8) 



which is a two-body correlator [36| . This remains finite in 
the TG limit g — > oo due to the simultaneous vanishing 
of the zero-distance density-density correlations and of 
the ID scattering length. Its expression in the TG limit 
can be found in |37j . 

We start again from Eq.Q to estimate the weight of 
the high-momentum tails. They are related to the short- 
distance non-analytic behaviour of the one-body density 
matrix induced by the presence of the delta-interactions. 
The task is greatly simplified by the fact that the only 
term contributing to the high-momentum tails in the se- 
ries for piB(x,y) is the j = 1 term, which upon short- 
distance expansion yields [43[ 



Pib I) ( a; )») 



\x - 



F(R), 



(9) 



where we introduced the center-of-mass coordinate R = 
(x + y)/2 and the two-body function F(R) is given by 

F(R) = n (R) MdRUu{R)\ 3 - \ £ f v u v (R)d R u* v (R)\ 2 

( 10 ) 

with n(R) = J2v fu\u u (R)\ 2 being the density pro- 
file. Using the property of the asymptotics of 
Fourier transforms, J Aze~ lk ^ z ~ z ^\z — z \ a ~ 1 F(z) — > 



C 



class 



N 2 



knT 



3 / 2 V huj 



(14) 



We thus find that a classical Tonks- Girardeau gas in har- 
monic confinement displays large-momentum tails with 
a weight which increases as the square-root of the tem- 
perature. On the other hand, at increasing temperature 
the wavevector fc starting from which the momentum 
distribution displays algebraic tails also increases with 
temperature, ie k oc 1/XdB hi the classical regime, with 
XdB = \/2irh 2 /rnksT '. We would also like to notice that 
the behaviour of the contact at large temperatures de- 
pends on the type of external confinement; for a box trap 
a similar calculation yields that the contact increases lin- 
early with temperature. 

The temperature dependence of the contact obtained 
from Eqs.(fl"0| and ([TTj) is illustrated in Fig 12 and com- 
pared to the results from numerical evaluation of the 
full momentum distribution and the high-temperature 
limit (HU). 

Virial approach for the contact at high temperature To 
gain further insight in the large-temperature behaviour 
of the contact, we derive it using a virial approach. 

The finite-temperature contact can be obtained from 
from Tan's sweep relation [a, |24| in its thermodynamic 
form 39] 



dan 



2m 



(15) 



where a\D = —2H 2 /mg and Q is the grandthermody- 
namic potential. Eq. (|15[) follows from the Hellmann- 
Feynman relation together with the thermodynamic 
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ksT /huj 



FIG. 2: (Color online) Tan's contact (in units of a\ a ) as a 
function of temperature feaT (in units of ftw) for a TG gas un- 
der harmonic confinement. The full expression from (|10p and 
(JTTJ (solid, magenta) is compared with the high-temperature 
limit (|14l) (dashed, red) and the data from the numerical cal- 
culation of the momentum distribution (crosses, blue). 



identity relating energy and grandthermodynamic poten- 
tial [39(. The advantage of this formulation is that it 
yields a virial expansion form for the contact. 

We start from the virial expansion of the grandther- 
modynamic potential ft — —ksTQi{z + b2Z 2 + b^z 3 + ....) 
in powers of the fugacity z — e 13 ^, with & 2 = tt- 



2 



Ql 

>2 = Qi J2u e ~^ e " and more generally Q n = Tre~ 
is the n-body cluster. The energies e]f entering in b-i 
correspond to the eigenvalues of the two-body problem 
in the relative-coordinate variable. Using the Tan's sweep 
relation, the above virial expansion turns then into the 
virial expansion for the contact, 



C 



2m 



dB 



■k B TQx{c2Z 2 + c 3 z 3 + ...) 



with adimensional coefficients given by 

db n 



d(a 1D /X d B) 



(16) 



(17) 



As a consequence of the scale invariance associated to the 
Tonks-Girardeau limit we can predict a universal feature 
of the coefficients c„ , namely we expect that they are con- 
stant (ie independent on temperature) since (oid/A^b) 
vanishes in the TG limit. 

The direct evaluation of C2 in the TG limit follows. In 
harmonic trap we have that e]f- = fkx)(v + 1/2) with v 
given by the solution of the transcendental equation 



r(-v/2) 



V2au 



r(-i//2 + l/2) a ho 



(18) 



The TG regime corresponds to am = 0, hence v = 2n+l 
with n = 0,l,2,3... and the derivative required in (|17|) is 



readily evaluated. We finally obtain 

2 3 / 2 /3fiwAds 

n=0 



C2 



r(n + 3/2) e _ J 9ft a) (2n+3/2) 



1.! 



(19) 



The sum is performed exactly to obtain 



C2 



i for fifa < 1. (20) 
v2 



We hence find a constant value of the two-body contact 
coefficient at large temperature, as expected from uni- 
versality considerations. Using the above result for c 2 , 
and the classical expression for Qi = ksT/hui and for 
the fugacity z = NUtt/hsT, we readily recover Eq. (|T4"]) . 

Conclusions. We have developed the formalism of the 
thermal Bose- Fermi mapping to calculate the momentum 
distribution of a Tonks-Girardeau gas under harmonic 
confinement at finite temperature. As a main feature, we 
have found that its high-momentum tails have a weight 
- the Tan's contact - which increases with temperature, 
and have linked this to the effect of reduced dimension- 
ality. Moreover, we have shown that this result can be 
understood as one of the signatures of universality in TG 
gases, associated to scale invariance in the infinite inter- 
action limit. Our results seem within experimental reach 
and open the way to exploring further aspects of 
universality in strongly interacting ID gases, both homo- 
geneous and under confinement. 
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SUPPLEMENTARY MATERIAL FOR: UNIVERSAL CONTACT FOR A TONKS-GIRARDEAU GAS AT 

FINITE TEMPERATURE 

Derivation of the fermionic expression for the bosonic one-body density matrix 

We derive here the expression (4) for the bosonic one-body density matrix at finite temperature in terms of fermionic 
quantities. As a first step, due to the effect of the mapping function in Eq.(2) of the main text each many-body integral 
in Eq.(3) of the main text has the form 

/ dx 2 , ...dx N ^N,a(x,x 2 ..,x N )^* N a (y, x 2 , --,x N ) = TT / da; l sign(a; - x t ) sign(y - Xi)f(x, y, x 2 , ...x N ) (21) 

where for short-hand notation we have set f(x,y,x 2 ...X]y) = vpjy a (x, x 2 ...XN)^f ! * a (y, x 2 ...xn)- We use that for each 
integration variable Jj dxisign(x — Xi) sign(y — xi) f = Jj dxi f — 2 sign(ir — y) J x dxi f and the binomial power series 
to rewrite Eq.(3) of the main text as 

/TV — l\ 

p 1B (x,y) = ^2Pn,uN^2 ( • ) (- 2 ) 1 ( s[ gn(x - y)) 3 / dx 2 ...dx j+1 

N,a j=0 \ 3 / 

x J dx 0+2 ...dx N ^x a (x, x 2 .., xjv)*^* Q (y, x 2 , .., x N ). 
Here, we recognize the fermionic j-body density matrices, which are defined as 

p jF (x 1 ..x j ;x' 1 ..x' j ) = y^Pjv,q . N ' ,. ( / dx j+1 ..Ax N ^^ a (x 1 ..x 3+1 ..x N )^^ l * a (x' 1 ..,x j+1 ..,x N ). (22) 

N,a ^ 3)- Jl 

Thus, we are led to the expression for the bosonic one-body density matrix as a sum of integrals of fermionic density 
matrices, 



~ r_ 2)T . ry 

Pib{x,v) = 22 — —( s ig n ( x ~y)y dx 2 ...dx j+ ip J+ i tF (x,x 2 ,...x J+1 ;y,x 2 ,...x j+ i). 

3=0 J ' Jx 



(23) 
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Further progress can be made by noticing that for noninteracting fermions the j-body density matrices factorize in 
terms of the corresponding one-body density matrices according to 

p jF (x 1 ,x 2 ,...x n ;x' 1 ,x' 2 ,...x' n ) = det[p 1F (xi,x' e )]i,e =hn , (24) 

where pip(x,y) — J2f=i fvj u vj{x)u%.{y) and /„ = \/[eP( £ «~^ + 1]. The final expression for the bosonic one-body 
density matrix (4) in the main text is thus obtained. 

Matrix functions B ViV . entering the determinantal expression of the one-body density matrix under 
harmonic confinement 

In the main text we have demonstrated how the initial multi-variable integration needed to compute the one-body 
density matrix as given by Eq. (4) of the main text can be reduced to the calculation of the minor determinant of 
the matrix B(a;, y). Here we give its full expression for the case of a TG gas under external harmonic confinement. 

As it is readily obtained from Eqs.(5) and (6) of the main text, the expression for the matrix B(x,y) is given by 

BviVj (x, y) = sign(z - y) / dwu Vi (w)ul.(w). (25) 

Jy 

For harmonic confinement, the calculation of the one-body density matrix can be further simplified by the fact 
that the elements of the matrix B ViV .{x,y) can be evaluated analytically. Using the power series expansion of 

the Hcrmitc polynomials H v {x) = v\ J2kio {-l) k (2x)"- 2k /k\/(v - 2k)\ and the result p m (x) = J x dwe~ w ' '' w m = 
x m+i e ~x p^i^ (m + 3)/2, x 2 ) where iF\(a,b,c) is the confluent hypergeometric function, we have 

l 1"«/2] to/2] (_nfc+k' 

>■"'■■ !/) = -H"^ E E fe! ^!( y - 2fe)!( y - 2feO! ^ + ^- 2(fc+fc ' )(:C ' y) - (26) 

I 



